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Abstract 

We study a stiff quasi-periodic orbit of the electromagnetic two-body problem of Dirac's elec- 
trodynamics of point charges. The delay equations of motion are expanded about circular orbits 
to obtain the variational equations up to nonlinear terms. The three-frequency orbit involves two 
harmonic modes of the variational dynamics with a period of the order of the time for light to 
travel the interparticle distance. In the atomic magnitude, these harmonic modes have a frequency 
that is fast compared with the circular rotation. The quasi-periodic orbit has three frequencies; 
the frequency of circular rotation {slow) and the two fast frequencies of two mutually orthogonal 
harmonic modes. Poynting's theorem gives a mechanism for a beat of the mutually orthogonal fast 
modes to cancel the radiation of the unperturbed circular motion by interference. The nonradia- 
tion condition for destructive interference is that the fast frequencies beat precisely at the circular 
frequency. The resonant orbits have magnitudes in qualitative and quantitative agreement with 
quantum electrodynamics (QED), as follows; (i) the orbital angular momenta are integer multiples 
of Planck's constant to a good approximation, (ii) the orbital frequencies agree with a correspond- 
ing emission line of QED within a few percent on average, (iii) the orbital frequencies are given by 
a difference of two linear eigenvalues, viz., the frequencies of the mutually orthogonal fast modes, 
and (iv) the angular momentum of gyration of the variational motion about each resonant orbit is 
of the order of Planck's constant. 
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I. INTRODUCTION 



We study a stiff quasi-circular orbit of the electromagnetic two-body problem of Dirac's 
electrodynamics with retarded-only fields[l], a dynamics with implicitly-defined delay. The 
motivation is to understand this complex dynamics described by the delay equations for 
particle separations in the atomic magnitude[2]. We give an economical method to derive 
the variational equations of dynamics about circular orbits up to nonlinear terms. The 
particular harmonic solutions of the variational dynamics, with a period of the order of the 
time for light to travel the interparticle distance, are henceforth called the ping-pong modes 
(PP). For atomic orbits, the frequency of the PP modes turns out to be much faster than the 
orbital frequency. These modes introduce a fast (stiff) timescale in the dynamics, and are 
physically important for the particles to have the option to avoid radiating energy away. The 
PP modes for vibration along the orbital plane turn out to have almost the same frequency 
of PP modes for vibrations normal to the orbital plane, a remarkable quasi- degeneracy 
that naturally produces beats with a frequency of the order of the orbital frequency. The 
quasi-periodic orbits of hydrogen have three frequencies; the slow frequency of the circular 
orbit plus the two fast frequencies corresponding to a planar PP mode and a perpendicular 
PP mode, as illustrated in Fig. 1. This special combination appears in a mechanism 
suggested by Poynting's theorem to cancel the radiation of the slow circular motion by 
interference with a beat oscillation of the two mutually orthogonal PP modes. We investigate 
the conditions for the nonlinear variational equations to accept such fast gyrating solutions 
about the circular orbit, i.e., the PP oscillations. After the fast dynamics is established, 
Poynting's theorem gives a necessary resonance condition to avoid radiative losses; viz., the 
two mutually orthogonal PP modes must beat at the orbital frequency. This resonance 
condition turns out to be satisfied precisely in the atomic magnitude. The fast gyration 
defines an angular momentum vector of the order of the orbital angular momentum of 
the unperturbed circular orbit. We stress that the point charges are not spinning about 
themselves, but rather gyrating about a guiding-center that is moving along a slow circular 
orbit, as illustrated in Fig. 1. The stiff three- frequency orbits share several magnitudes 
with those of the hydrogen atom of quantum electrodynamics (QED) [3], with reasonable 
precision and qualitative detail. The circular frequency of a resonant orbit agrees with the 
corresponding line of QED within a few percent average deviation. There is also a large body 
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of qualitative agreement with QED; (i) the resonant orbits have orbital angular momenta 
that are approximate integer multiples of a basic angular momentum. This basic angular 
momentum agrees well with Planck's constant, (ii) the angular momentum of gyration of 
the variational motion about the circular orbit is of the order of Planck's constant. This 
angular momentum of gyration is a vector that rotates at the orbital frequency, and (iii) the 
emitted frequency is given by a difference of two linear eigenvalues, i.e., the frequencies of 
the PP modes, analogously to the Rydberg-Ritz principle of QED. 

The equations of motion of Dirac's electrodynamics of point charges[l] are briefly dis- 
cussed in Appendix A. After Dirac's 1938 work[l], an early study of Ehezer [4-7] revealed 
the surprising result that an electron moving in an attractive Coulomb field can never fall 
into the center of force by radiating energy (henceforth called Eliezcr's theorem). It was 
subsequently found that only scattering states are possible in any tridimensional motion 
with self- interaction in a Coulomb field [6, 7]. Ehezer's theorem strongly suggests that a 
finite mass for the proton is essential for a physically meaningful dynamics in the electromag- 
netic two-body problem. When the proton has an infinite mass, there is an inertial frame 
where it rests at all times, and in this frame the protonic field on the electron is simply a 
Coulomb field, i.e., the dynamics of Eliezcr's theorem. On the other hand, if the proton 
has a finite mass, such inertial frame does not exist and the equations of motion involve 
delay, because of the finite speed of light. A finite mass for the proton is what brings delay 
into the electromagnetic two-body dynamics, with its associated ping-pong phenomenon. 
The infinite- mass limit is a singular limit, because the equations of motion pass from delay 
equations to ordinary differential equations. This work is an attempt to put together what 
is lost in this singular limit where the PP modes disappear. Wc stress that the circular orbit 
is not an exact solution of the full equations of motion, so that we are not doing Lyapunov 
stability, but rather constructing particular solutions of the nonlinear variational equations, 
i.e., the three-frequency orbit. 

The road map for this paper is as follows. In Appendix A we review the electrodynamics 
of point charges in a generahzed setting that includes Dirac's theory as a special case, and 
give the equations of motion of point charges in an intuitive form. A non-specialist reader 
should start reading the paper from Appendix A. In Section II we define the PP modes and 
the quantities of the circular orbit, to be used as an approximate solution in Sections III, IV 
and V. In Section III we outline our economical method to derive the variational equations. 
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a method that expands the imphcit hght-cone condition and uses the action formahsm. In 
Section IV we derive the hnearized variational equations and study the PP modes of tangent 
dynamics for vibrations along the orbital plane. This derivation is laborious and makes full 
use of our economical method plus the use of a symbolic manipulations software. We derive 
the linearized variational equations in the generalized electromagnetic setting of Appendix 
A, to compare with previously known results, but we stress that in Section V we use only 
Dirac's electrodynamics, the physically interesting special case. In Section V we give an 
application of Dirac's electrodynamics to atomic physics, by discussing the existence of a 
three-frequency orbit involving mutually orthogonal PP oscillations of finite amplitude, i.e., 
a particular solution of the variational equations. Wc investigate the mechanism to cancel 
the radiation of the circular orbit by interference with a beat of two mutually orthogonal PP 
modes, a mechanism that starts to operate immediately after the PP dynamics is established. 
We discuss the necessary resonance condition of Poynting's theorem to avoid radiation, i.e., 
that the PP modes beat at the orbital frequency. We study the resonant orbits that are 
stabilized by this mechanism and compare their magnitudes with the magnitudes of QED. 
In Section V we also give a second derivation of Poynting's resonance condition. This 
derivation averages the angular momentum of gyration over the fast timescale, yielding a 
vector rotating at a slow frequency. The rotating angular momentum produces a gyroscopic 
torque on the slow dynamics, and introduces the same resonance condition of Poynting's 
theorem. In Appendix B we derive the tangent dynamics for oscillations perpendicular 
to the orbital plane, analogously to what is done in Section IV for the planar variational 
equations. These two derivations can be given separately up to the linear order, but the z 
and xy variational equations are otherwise coupled at higher order. The existence of fast 
harmonic solutions of finite amplitude is discussed in Appendix C. In Appendix C we also 
give a third derivation of the resonance condition, a derivation based on the detailed balance 
of the guiding-center dynamics. Last, in Section VI we put the conclusions and discussion. 



II. THE CIRCULAR ORBIT 



Our perturbation scheme takes the circular orbit as a first approximation. We use the 
index i = 1 to indicate quantities of the electron and i = 2 for the proton, with masses m\ 
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and 777.2, respectively. We henceforth use a unit system where the speed of hght is c = 1, 
the electronic charge is ei = —62 = — 1 and the mass of the electron is 777i = 1. In our unit 
system, the mass of the proton is given by 7772 — 1824, approximately. The circular orbit is 
illustrated in Fig. 2; the particles move in concentric circles and in diametral opposition at 
the same time of the inertial frame. The details of the familiar Coulombian circular orbit will 
be given now. The constant angular velocity is indicated by O, the interparticle distance in 
light-cone is r;, and the angle that one particle turns while the light emanating from the other 
particle reaches it is ^ = Jlr;,. The delay angle Q is the natural independent parameter of 
the circular orbit, which turns out to be small for orbits of the atomic magnitude, Q < 10~^. 
For small the interparticle distance in light-cone, r^,, is 0{(P') close to the interparticle 
distance at the same time of the inertial frame, Tq. Because of this Oip^) approximation, 
the familiar Coulombian formulas with replaced by ro yield the leading order formulas 
in powers of Q. For example, the orbital frequency, Kepler's law, is given to leading-order in 



where // = 77717772/ (7?7i -|-7?72) is the reduced mass (here and henceforth). In our unit system, 
// = (1824/1825) ~ 1 for hydrogen. The interparticle distance in light-cone, r^, is constant 
along the circular orbit and given, to leading order, by 



Using the radial equation of motion for the Coulombian orbit, one finds that the angular 
momentum of the circular orbit is given, to leading order in ^, by 



The units of 1^ in Eq. (3) are e^/c, just that wc arc using a unit system where = c = 1. 
For orbits of the atomic magnitude, 1^ as defined by Eq. (3) is of the order of one over the 
fine-structure constant, aT^ — 137.036. Each particle travels a circular orbit with radius 
and scalar velocity defined by 



^by 



(1) 




(3) 



(4) 



T2 = &2r6. 
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and 



= nn = 9bi, (5) 
V2 = fir2 = 6b2, 

for the electron and for the proton, respectively. For consistency with the definition of rj, 
as the exact interparticle distance in light-cone, the definition of 61 and 62 niust include an 
0(9^) term; 61 = (1 + g6^)m2/ (mi + 1112) and 62 = (1 + (^1 + ^2)- Evaluating the 

interparticle distance in light-cone and equating it to yields 

rl^rl{bl + bl + 2b,b2Cos{e)). (6) 

Formula (6) becomes Eq. (3.1) of Ref.[8] after use of Eq. (5). Expanding Eq. (6) for small 
6 yields g 0.5/x/(mi + 1712). The ratio (61/62) = {m2/mi) is the Coulombian ratio and 
the radii defined by Eq.(4) arc 0(9'^) near the Coulombian radii, and we henceforth call this 
the Coulombian circular orbit. As discussed in Appendix C, for Dirac's electrodynamics 
the circular orbit is only an approximate solution. Since we are doing perturbation, it 
suffices to use the above defined Coulombian orbit as an approximate solution, and let the 
perturbation scheme take care of the correction. In the action-at-a-distance theory [8, 9], it 
turns out that the circular orbit is an exact orbit of the two-body problem. The action-at- 
distancc electrodynamics is a special setting described by the equations of Appendix A with 
r = —1/2, and is used here only to crosscheck the method of Section III. 

The intuitive picture of a PP oscillation is a ping-pong game, i.e., the particles throwing 
a ball back and forth at the finite speed c — 1, as a means to communicate changes in 
position. The PP modes have a period of the order of r^/c, where Vb is given by Eq. (2). 
The order of magnitude of the PP frequency is 

wpp — T — 2TTa9^, (7) 

n/c 

where we used that c = 1 in our unit system. To compare the PP frequency with the 
Coulombian orbital frequency (1), it is useful to express the PP frequency (7) as 

wpp ~ 27ri^9^ = (8) 
9 

where |A| is a number of the order of 27r. For atomic orbits, 9 ~ 10~^, the PP frequency 
Wpp is larger than Q by three orders of magnitude. We henceforth define the generafized 
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frequency of a ping-pong normal mode by Xft/9, so that a purely imaginary A defines a 
harmonic oscillation, and the eventual real part of A defines a damping or a runaway. For 
small 9, the bouncing-time for light to travel back and forth the interparticle distance is 
approximately 2ri,/c, so that the phase-shift of the fast oscillation during the bouncing-time 
is A$ = Im(^^). This phase-shift evaluates to A$ = 2 lm(A) by use of c = 1 and f2rf, = 9, 
and we shall see in Section V and in Appendix C that A must be an integer multiple of 27r, 
i.e., X — irq. The integer number q is henceforth called the ping-pong index. 



III. VARIATIONAL EXPANSION ABOUT A CIRCULAR ORBIT 



The variational equations are obtained by substituting a circular orbit plus a perturbation 
into the equations of motion of Appendix A, e.g., Eq. (61), and expanding in powers of the 
perturbation. Since these are complex equations, even after using a convenient coordinate 
system this derivation is long. Our method uses the fact that the equations of motion are 
almost the Euler-Lagrange equations of a suitable Lagrangian, with the addition of the 
self-interaction force. After the variational equations of the Lagrangian sector are derived, 
we add in the variation of the terms due to the self-interaction. Readers should consult 
Appendix A for a brief review of the equations of motion of point charges and definition of 
terms such as the Lorentz-force, the Lorentz-Dirac self-force and the action integral. In this 
Section we outline the procedure of substituting a circular orbit plus a perturbation into 
action (68) of Appendix A and expanding the action up to a desired order. Minimization 
of this truncated action plus the expansion of the kinetic action yields the Lorentz sector 
of the variational equations. Since we perform the algebraic computations with a symbolic 
manipulations software, it is equally easy to derive the variational equations in a generalized 
electromagnetic setting that contains Dirac's electrodynamics as a special explained 
in Appendix A. The calculation in the generalized setting has some interest of its own, and 
provides a useful cross-check with other known calculations [2, 11]. The generalized setting 
of Appendix A has an arbitrary constant T in the Green's function, and we stress that F will 
be set to zero in our application of Section V, i.e., we study the hydrogen orbit in Dirac's 
electrodynamics with retarded-only fields [1], a fundamental physical theory. A reader not 
interested in this generality can assume F = throughout the whole paper. 
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The variational equations for planar perturbations of circular orbits decuple from the 
equation for transverse perturbations up to the linear order. It is convenient to write these 
planar variational equations using complex gyroscopic coordinates rotating at the frequency 
fl of the unperturbed circular orbit. The coordinates {xk, Uk) of each particle are defined by 
two complex gyroscopic coordinates r]^ and ^j. as 

Uj = Xk + iyk = nexp{int)[dk + 2r]^], (9) 
u* = Xk- iyk = rbexp{-int)[dk + 2^^;.], 

where A; = 1 for electron and A; = 2 for proton and the real quantities di = bi and d2 = —1)2 
are defined in Eqs. (4). Because Xk and yk are real, we must have rj^. = ^l, but to obtain the 
variational equations it suffices to treat rj^^ and ^j. as independent variables in Lagrangian 
(69). Two quantities appear so often in the calculations that it is useful to name them; (i) 
The numerator of Lagrangian (69), evaluated along the circular orbit, is constant and given 

by 

C= 1 + 6162^^ cos(^), (10) 

having the same value for both retarded and advanced interactions, and (ii) The denominator 
of Lagrangian (69), evaluated along the circular orbit and divided by r^, is constant and 
defined by 

5 = 1 + 6iMsin(^), (11) 

having also the same value for both retarded and advanced interactions. For the stiff hmit of 
Sections IV and V we shall ignore the 0{9'^) corrections and set C = 5" = 1. Here we derive 
the electron's equations of motion only (j = 1 in Eq. (9)). The equation for the proton is 
completely symmetric and can be obtained by interchanging indices 1 and 2. Notice that 
the exchange operation on the d's is di <^=^ — ^2- This is because di = bi is defined positive 
while d2 — —62 is defined negative, so that at the same time the particles are in diametral 
opposition along the unperturbed circular orbit, as illustrated in Fig. 2. 
The velocity of the electron at its time ti is calculated with Eq. (9) as 



ill = vix + iviy = 9 exp{iD,ti)[idi + 2i{'r]i — irji], 
ill = '^is ~ ^f^ij/ = ^exp(— iQii)[— idi — 2i{^i + 



(12) 



The velocity of the proton at time t2 can be obtained by interchanging indices 1 and 2 in 
Eq. (12), as explained above. Using the c — ±1 convention explained above Eq. (70) of 
Appendix A, the quantities of particle 2 appearing in Lagrangian (69) are evaluated at a 
time t2 = t2c, in light-cone with the position of particle 1 at time ti. The implicit light-cone 
condition must be expanded and solved by iteration, and for that it is convenient to define 
in each case the excess-lag function (p^ by 

t2c^tl + - + ^. (13) 

If the perturbation is zero, then (p^ — (p_ — and the dynamics is along the original 
circular orbit, where the light-cone lag is the constant for the advanced case and — 
for the retarded case. We henceforth measure the evolution with the scaled-time parameter 
r = Qti. The implicit definition of y?^ by the light-cone condition involves the position of 
particle 2 at either the advanced time t2+ or the retarded time ^2- , as defined by Eq. (9), 

U2{t + c9 + ifj = nexp{int2c)[d2 + 27]^{t + 00 + (^J], (14) 
uI{t + c9 + cpj = r6exp(-if]t2c)[c?2 + 2{2(''" + + V'c)]) 
as well as the velocity of particle 2 at the advanced/retarded position; 

U2{t + c9 + (/? J = i9 ex.p{iQt2c) [^2 + (% - + '2Vc{V2 - V2)] > (15) 
u*{t + c9 + (p,) = -i9 exp{-int2c) [d2 + {C2 + + 2<^c(^2 + ^2)] • 

To obtain the linearized variational equations we expand the equations of motion to first 
order in 77 and ^^j., which are the Euler-Lagrange equations of the quadratic expansion of 
action (68) in rj^ and We must therefore carry all expansions up to the second order in 
the Tj^ coordinates. For example, the position of particle 2 is determined up to the second 
order by expanding the arguments of and ^2 of Eq. (14) in a Taylor series about the 
unperturbed light-cone for one order only as 

U2{t + c9 + (p^) ~ n ex.p{iQt2c){d2 + 2[r]2{T + c9) + ipj^^i'r + c9)\}, (16) 
uI{t + c9 + if^) ~ nexp{-i^lt2c){d2 + 2[e2(T + c9) + tp^^^^r + c9)]}. 

In the following we find that ip^ is linear in 77 and ^ to leading order, so that the next 
term in expansion (16) would be a third order term, not needed for the linear variational 
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equations. We henceforth indicate ^77 quantities of particle 2 evaluated at the unperturbed 
light-cone by a subindex c, i.e., 772c = V2{'^ ''<' ^) = C2(''' + c^)- Otherwise when 

the subindex c appears, it indicates evaluation on the exact light-cone. Notice that the 
small parameter of expansion is the size of rj and ^, and we henceforth expand any quantity 
evaluated on the light-cone in a Taylor series about the unperturbed light-cone, up to the 
order needed. For quasi-circular orbits this method yields variational equations with a fixed 
delay. We stress that one should never expand the arguments in powers of 9; such an 
expansion yields ordinary differential equations and looses the PP modes! The perturbed 
light-cone is expressed implicitly by the distance from the advanced/retarded position of 
particle 2 to the present position of particle 1, described in gyroscopic coordinates by the 
modulus of the complex number 

U = Ui(T)-U2(T + c6 + ip^), (17) 

where c = 1 describes the advanced light-cone and c = — 1 describes the retarded light-cone. 
Using Eq. (16) to calculate u up to the second order yields 

u = rbex.p(int2c){D* + 2[exp(-ic^ - iip^)r]^{r) - r^^c - </'c^2c]}, (18) 
u* = rbexp{-int2c){Dc + 2[exp(ic^ + - ^2c - ^0^20]}^ 

where we defined the following complex function of (p^ 

Dc = b2 + biexp{ic9 + iLpJ. (19) 

At (p^ — (the unperturbed circular orbit), Dc has a unitary modulus, expressing the 
unperturbed light-cone condition (6). The right-hand side of Eq. (18) is a quadratic form 
times exp{iflt2c), and in the action it appears multiplied by a counter- rotating term, i.e., 
a quadratic form times exp{—iQt2c) , so that the product is independent of t2c- Because of 
this rotational symmetry of action (68), the following quadratic Gauge simplification can 
be applied directly to any quadratic rotating form. One can integrate by parts a quadratic 
term of the quantity, e.g., equation (18), and disregard the boundary term. This Gauge 
simplification yields a correct action up to the second order. For example, integrating by 
parts the quadratic terms in r/g^ and ^2c the right-hand side of Eq. (18) and disregarding 
the quadratic Gauge yields 

u = r^eMi^hcHD: + 2[(1 - iif,) exp(-ic^)77i(T) - (1 - ^M}, (20) 
u* = ri,eM-i^t2c){Dc + 2[(1 + i(^J exp(ic^)ei(T) - (1 - ip,)C2c]}, 
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where we have also expanded exp{i(f^) up to the hnear order in (p^, enough to give the correct 
quadratic action. Analogously, the velocity of particle 2 has the following expansion up to 
a quadratic Gauge 

U2{t + C0 + cp,) ~ ie&q){i^hc)[d2 + 2(%e - i%c)(l - <^c)], (21) 
u;{r + ce + ip,) ~ -^^exp(-^l]t2c)[(i2 + 2(^20 + 4c)(l-<^c)]- 

Using the above quantities, the numerator of Lagrangian (69) can be calculated as 

/l2 = (1 - Vi • V2c) = (1 - ^UlU*2 - ^U*i'^2), (22) 

and the denominator of Lagrangian (69) can be calculated as 

/i4 = ri2c(H ) + + ^ + (2^) 

In Eq. (23) we have introduced the scaled delay function (p by 

= ce<i>, (24) 

where the subindex under is omitted for simplicity of notation. To relate (j) to the ^rj 
perturbations we expand the implicit light-cone condition of the perturbed orbit up to the 

quadratic order 

uu*^\x2{ti + - + ip,)-xi{ti)\^ = c\t2c-h)'' = {n + ^y. (25) 
c lie 

The light-cone condition (25) is most simply expressed in terms of the scaled (f) defined in 
Eq. (24). This expansion of Eq. (25) up to the second order yields 

+2[{b2^^ - hiri2) exp(ic6') + (627^1 - 61^2) exp(-ic6')] 

+4[Ci?7i - Ci?72 exp(ic6') - ^2^1 exp(-ic6')] 

-2^[62(C2 + V2) + fciC2exp(-ice) + 6i972exp(ice)]. (26) 

The solution to Eq. (26) up to the first order in the ^77 coordinates is 
50(1) = [(61^1 - 62^2) + (^'1^1 - ^2%) + (^^2^1 - ^'1^2) exp(ic^) + (627?! - 617^2) exp(-ic^)]. (27) 
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Last, as a check for the above calculations, in the following we derive the equations of 
motion for the circular orbit of the action-at-a-distance theory [8]. The Lagrangian of action 
(68) for r — —1/2 is £ = (L_,_ + L_)/2, and its expansion up to the first order is 

2 = 1- {[e\S - 1)S + C^]{b^ + 62COs(e)) + S{esm{e) - e^cos{e))b2} ^'^'^J'-\ (28) 

where the tilde indicates that £ was scaled by its value along the unperturbed circular orbit, 
£ = rjjS £/C. Scaling the kinetic Lagrangian with the same factor and expanding to first 
order yields 

Ti - + ^ iVi + ^i), (29) 



where 7^ = l/^/l — vf and Vi is given by Eq. (5). The effective Lagrangian for particle 1 
up to the linear order is 

L« ^f, + £. (30) 

Lagrangian (30) is a linear functional oi ^i, independent of ^1, so that the Euler-Lagrange 
equation for is simply — 0, i.e., 

mibin^^e^s^ = [c^ + e^s{s - 1)] (61 + 62 cos(^)) + s{e sm{e) - cos{e))b2. (31) 

This is Eq. 3.2 of Ref. [8], and the equation for rji is the same condition by symmetry (this 
is actually the reason why the circular orbit is a solution[13]). The equation for particle 2 
is obtained by interchanging indices 1 and 2 in Eq. (31), as is Eq. 3.3 of Ref. [8]. In the 
next section we expand the action to second order, to determine the linearized variational 
equations. 

IV. LINEARIZED VARIATIONAL EQUATIONS 

In this Section we obtain the linear-order terms of the variational equations. For this we 
carry the expansion of Section III to the quadratic order. The next term of expansion (29) 
of the kinetic Lagrangian of particle 1, calculated with Eq. (12), is 

Ti = To 1 - Pir = ^ (^1 + Ci) (32) 



2C 
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This kinetic Lagrangian of particle 1 has a quadratic form defined by two coefficients 



~ rb9^ma^Sbi 

Ti = ^ iVi + ^i) (33) 

where Mi = {1 + ^l)mi^irbO^ S / C and Gi = (7^ — l)mi^i9^S/C. We also need the solution 
of Eq. (26) to second order, = ^(i)+ ^(2) > where ^^^^ is given by Eq. (27) and 0^) is 
calculated by iteration to be 



5'0(2) = viihCi - ibir]2 - 617)2) exp(ic6l) + </'(i6iC2 - ^^2771 - &1C2) exp(-ic6l) (34) 

1 

+2[Ci^i - Ci^2exp(ic^) - C2^iexp(-ic^)]. 



Next we expand the numerator and the denominator of Lagrangian (69), Eqs. (22) and (23), 
up to the quadratic order. We also need the following quantities expanded up to second 
order; (i) the particle separation in light-cone, Eq. (18), (ii) the velocity of particle 1, Eq. 
(12), and (iii) the velocity of particle 2 at time Eq. (21). This quadratic form can be 
greatly simplified by adding suitable Gauge terms to it (essentially integration by parts), 
which yields the following normalized Lagrangian for particle 1 

Lc = -^i?2i(r?iei - eir/i) + t/iiCir;! + \Ni,il + (35) 

+^0^1^20 + R*M2c + Pciiri2c + P*cVli2c + 

Ac • A* • 

Y(Ci^2c - ^2cCl) + y(^iCi - C2'7l) + 

Tr rp:)f 771 771* 

In Eq. (35) the coordinates of particle 2 appear evaluated in either the retarded or the 
advanced unperturbed light-cone, as indicated by the subindex c. The coefficient of each 
normal-form binary is a function of ?7t.i, m2 and 9, and is obtained by a Gauge- invariant 
combination of derivatives, e.g., the coefficient i?2i is the same in both retarded and advanced 
interactions and given by 

13 



The coefficients were evaluated witfi a symbolic manipulations software, starting directly 
from the effective Lagrangian and taking the necessary derivatives (in this way we avoid 
mistakes). The explicit functional dependences are not given here for brevity, but were 
checked at the various limits; Biot-Savart field. Coulomb interaction, Darwin interaction, 
etc. The effective Lagrangian for particle 1 is composed of the partial Lagrangians (69), 
evaluated at the advanced and the retarded light-cones, as in Eq. (68), plus the kinetic 
Lagrangian, 

L^y^ = ri-rL+ + (i + r)L_. (37) 

The linearized Euler-Lagrange equation of Lagrangian (37) respect to is a linear function 
of coordinates Ci ^ Vi ^ C2+ ) V2+ > C2- V2-i ^ ^^^11 as of their first and second derivatives, 
i.e., 

hdCi, Vi, C2+, V2+, C2-, V2-) = (38) 
-[{Nn + e^G)C, + e^Gl,] + [Mi(7)i + 2ir), - rj,) - UnVi - iB2ir]i] 

-(i?;e2+ + ^1^2-) - (^^+^2+ + YM-) + {T+^2+ + TJ^2-) 
-{P+V2+ + P-V2-) - {KV2+ + A-'72-) + {E+V2+ + E_fj2_). 

The Euler-Lagrange equation of Lagrangian (37) with respect to rji, ^2 ^^'^ V2 yields three 
more linear equations, which together with Eq. (38) form a system of four linear delay 
equations. 

In the following we explain how to include self- interaction into Eq. (38). The Euler- 
Lagrange equation of the kinetic energy (33) with respect to can be expressed as 

(1 + ^i)maAS/C)rln\ri, + 2ifi, - rj,) (39) 

= rl{S/C) — {mi-f^rbUi), 

where we replaced 9^ by r^Jl^ into the definition of Mi given below Eq. (33). On the second 
line of Eq. (39) wc recognize the variation of the complex momentum, mi^i^ryiii = pix + ipiy, 
multiplied by the factor rl{S/C). The variation of momentum is the force, so that to account 
for self-interaction we add to Eq.(38) the x-component of force (71) multiplied by the factor 
rl{S/C), plus i times the y-component of force (71) multiplied by rl{S/C), i.e.. 
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where the dots represent derivative respect to scaled time t = D,t. Using Eq. (40) with ui 
given by Eq. (9), and expanding up to the hnear order, yields the order-zero offensive force 
(94) plus the linearization of the self-force in gyroscopic coordinates. The offensive force (94) 
is the nonhomogeneous term of the variational equation, and shall be dealt with in Appendix 
C. In this Section we discard it and keep only the linear part of the variational equations. 
Notice that Eq. (39) came out naturally in the form of variation of momentum multiplied 
by rl- This instructive normalization suggests that we scale the equations of motion with 
the size of the unperturbed Coulombian attraction, l/rf. This normalization is useful when 
discussing orders of magnitude, and it is used in Appendix C to discuss estimates with an 
intuition about physical orders of magnitude. 

The linear variational equations for rji, ^2 V2 form a set of four linear delay 
equations, a system that can be solved in general by Laplace transform [14]. In the following 
we focus on the planar normal-mode solutions of this linear system, with a fast ping-pong 
frequency Xxy^/0, a definition motivated by Eq. (8). The complex number X^y is so 
far arbitrary, but a harmonic oscillation is defined by an imaginary X^y. In the following 
we substitute = Aexp{Xxyflt/9), rji = B exp{Xxyflt / 9) , = C exp{Xxyflt / 9) and 7^2 = 
D exp{Xxy^t/9) into the linearized equations and assume 9 small and \Xxy\ of order-one or 
larger, as discussed below Eq. (8). This yields four homogeneous linear equations for A, B, 
C and D, and a nontrivial solution exists only if the determinant vanishes. Using a symbolic 
manipulations software this determinant evaluates to 



2{l + 2T)9^Xxy , (l + 2r)^% 2 /X 

i4 



1 _ zi^_^^rL^ + ^ ^ - ^^(1 + 2Tf9'Xl + ... (41) 



+^(1 + ^ + + 2r) sinh(2A.,) - 2(1 + 2r + cosh\Xxy)] 

xy xy 

-2^(^ + J-)[2(l + 2r) cosh2(A,,) - (1 + 2r + 2T^) sinh(2A,,)] = 0, 

where M = mi + 1712- In our unit system, the total mass of hydrogen is M = 1825. For 
Dirac's retarded-only electrodynamics, Eq. (41) with F = 0, we obtain the following planar 
normal-mode condition for ping-pong modes (the stiff-limit) 

2 7 10 5 , , 

(1 + 3- + + - + • • • ) (^) exp(-2A.,) (42) 

Axy Axy Ky ^xy 

2 2 1 4 2 
= ^- Ky + -^9 X^y^ , 
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Comparing Eq. (42) to Eq. (90) of Appendix B we find that the quasi-degeneracy phe- 
nomenon exists only in three cases; (i) F = 0, Dirac's theory with retarded-only interac- 
tions, (ii) r = — 1, a non-physical advanced-only interactions theory, and (iii) F = —1/2 
, the dissipative Fokker theory of Rcf. [2] and the action-at-a-distance electrodynamics[ll]. 
This discriminating degeneracy is an interesting curiosity, and in this paper we disregard the 
two other dynamics that exhibit the quasi-degeneracy phenomenon, (ii) and (iii). In the 
next Section we study a three-frequency orbit of the hydrogen atom in Dirac's electrody- 
namics with retarded-only fields (F = ), the physically sound choice to describe hydrogen 
in nature. 



V. THE THREE-FREQUENCY ORBIT 



As discussed in Section IV, in Dirac's electrodynamics there is a remarkable quasi- 

dcgcneracy of the perpendicular and the planar tangent dynamics. In the large-|A| limit, 
both Eq. (42) and Eq. (90) of Appendix B reduce to 

(^)exp(-2A) = l, (43) 

henceforth called the degenerate stiff-limit. The value of (fi/M) for hydrogen in Eq. (43) 
is a small factor of about (1/1825). Equation (43) requires that A has a negative real part 
given by Re(A) = —a = — \n{.^'^). For the atomic magnitude, ~ 137.036, the value 
of a is about a ~ 15.0. On the other hand, the imaginary part of A can be an arbitrary 
multiple of tt, so that the general solution of Eq. (43) is 

-a + nqi, (44) 

where i = and the integrer q is the ping-pong index. Notice that the real part of A 

is always negative, so that the tangent dynamics about the circular orbit is stable in the 
stiff- limit. The unfolding of the degeneracy comes with the terms of order 1/A^ and 6^X^, 
as found in Eq. (42) and Eq. (90) of Appendix B. The exact roots of Eqs. (42) and (90) 
near the limiting root (44) are defined, respectively, by 

Kyid) = -cr^y + Trqi + iei, (45) 
Xz{9) = —az + nqi + i€2, 
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where ei(^) and €2(0) are real 0{9) numbers. The three-frequency orbit is formed from an 
initial circular orbit as follows; (i) Initially, the self-force (94) dissipates energy, essentially 
the radiation of the circular orbit. Some of the radiated energy is absorbed directly by the 
PP oscillations. The slow guiding-center circular motion may also loose radius to account for 
some of the energy loss, spiralling in for a small number of turns, (ii) After the PP modes 
absorb enough energy, their amplitudes grow and the PP oscillations become neutrally 
stable, which is illustrated in Fig. 1. As discussed in Appendix C, at a finite PP amphtude, 
a harmonic solution to the nonlinear variational equations exists, i.e., the a's of Eq. (45) 
vanish, so that the A's become purely imaginary. This balancing is achieved at relatively 
small amplitudes, i.e., near the circular orbit, and (ill) The radiation of the PP modes starts 
to interfere with the orbital radiation if the orbit is a resonant one. In the following we 
discuss a surprisingly simple stabilization mechanism that operates after the fast harmonic 
oscillations are established near the guiding-center orbit. 

We henceforth assume that the three-frequency orbit of Fig. 1 is defined by the following 
combination of a planar PP mode and a perpendicular PP mode 

Xk + iyk = rbexp{int)[dk + 2Ck], (46) 

Xk - iyk = rbCxp{-iQt)[dl + 2^^], 

^k = 'l^b^k, 

with 

= R'','exp[{7rqt + te1)nt/e], (47) 
Zk = Re{Rlex-p[{7rqi + ie^2)^t/e]}, 

where R^^ and are respectively the amplitude of the planar PP mode and the amplitude 
of the perpendicular PP mode { k = 1 for the electron and k = 2 for the proton). Notice 
that we chose the same q for the two perpendicular oscillations of Eq. (47) , so that q cancels 
and the difference between the frequencies of Eq. (47) defines a beat in a slow timescale. 
The far-electric field of the electron, at a far distance Too, is obtained by exchanging indices 
and disregarding the first term on the right-hand side of Eq. (63) of Appendix A, 

n X (n X ai_) n x (v,_ x ai-) 
(1 - n • Vi_)'^roo (1 - n • Vi_)'^roo 
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where we included a minus-one factor to account for the negative electronic charge. In Eq. 
(48) , the unit vector n points from the retarded position of the electron to the observation 
point at infinity, while Vi_ and ai_ are respectively the Cartesian velocity and acceleration 
of the electron (particle 1). The Poynting flux is proportional to |Eip, so that the common 
denominator (1 — n ■ Vi_)^on the right-hand side of Eq. (48) can be factored off |Ei|, and 
is henceforth ignored. For the radiation of a circular orbit of atomic magnitude, the second 
term on the right-hand side of Eq. (48) is small and can be disregarded. This is the dipole 
approximation of atomic physics, valid when the interparticle distance is smaller than the 
radiated wavelength. This approximation is not valid for the three-frequency orbit of Fig.l 
because of the short wavelength of the PP modes. For the far-field of the three-frequency 
orbit of Fig. 1, the second term on the right-hand side of Eq. (48) must be kept, and is 
henceforth called the spin-radiation field. Because of this quadratic spin-radiation field, the 
vector product of the mutually orthogonal PP modes (46) beats at a slow frequency (cosine 
times cosine averages to cosine of the difference). The spin-radiation field of the electron, 
Epp, averaged over the fast timescale and assuming \Rf \ ~ ~ in Eq. (47), has the 
following amplitude 

IEppI = cos([l + fci^]l^t). (49) 

The oscillation of Eq. (49) is determined by the frequency difference (beat) of the mutually 
orthogonal PP modes of orbit (47) , 

Awpp = [1 + ^ ^ ^ '' ]n. (50) 

Because of the quasi-degeneracy property, 62 and are 0{d) quantities, so that Awpp in 
Eq. (50) defines a beat of the order of the guiding-center frequency Q. The first term of the 
right-hand side of Eq. (48), averaged over the fast timescale, yields the usual dipole field 
of the unperturbed circular motion, with magnitude 

|Ed| = — cos(nt). (51) 

Poynting's theorem gives a necessary condition for Epp to destruct E^ by interference, viz., 
both fields must oscillate at the same frequency. This yields the resonance condition 

n = Awpp. (52) 

Notice that at resonance, the three-frequency orbit defined by Eqs. (46) and (47) has only 
two independent frequencies [16], because of relation (52) between its frequencies. In the 
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gyroscopic coordinates (47), the resonant orbit has only one frequency, and in Appendix C 
such resonant orbit is called a harmonic solution. Use of Eq. (52) with Eq. (50) yields 

- e? = 0. (53) 

In the following we approximate el and £3 respectively by the ei and 62 of the linear problem, 
as defined by Eq. (45), assuming that essentially the real parts of the hnear eigenvalues (45) 
are modified by the finite- amplitude correction (93), while the imaginary parts acquire only 
a small correction. According to this approximation, we expand the small correction in 

powers of 6 and take the first term 

62-61 ^e^-e?+p^ + --- . (54) 

This approximation is justified in Appendix C and the linear coefficient p should depend on 
the radius of gyration, < l/nq. Condition (53) together with Eqs. (50) and (1) gives Q 
as 

n = iie\e2-ei) + {i-p)e^ + --- , (55) 

a Rydberg-Ritz-like formula. Resonant frequency Q is expressed by an eigenvalue difference 
in another way; the exact difference between the frequencies of the orthogonal PP modes, 
Eq. (52). The PP frequencies should be rigorous linear eigenvalues of the tangent dynamics 
of the three- frequency orbit, but here we consider only the neighborhood of the circular 
orbit, and keep to approximation (55). 

We henceforth set p = 1 into Eq. (54) as a qualitative approximation, yielding 

6i - 62 + ^ = 0, (56) 

which has the solutions listed in Table 1. The root-searching problem of Eq. (56) is well 
posed and for each integer q there is a single solution 9 to Eqs. (56), (42) and (90), i.e., 9 is 
quantized by the integer q. According to QED, circular Bohr orbits have maximal angular 
momenta and a radiative selection rule ( Al — ±1) restricts the decay from level A; -|- 1 
to level k only, i.e., circular orbits emit the first line of each spectroscopic series (Lyman, 
Balmer, Ritz-Paschen, Brackett, etc.), the fourth column of Table 1. We have solved Eq. 
(56) together with Eqs. (42) and (90) using a Newton method in the complex A-planc. The 
numerically calculated angular momenta 9~^ are given in Table 1, along with the orbital 



19 



frequency in atomic units (137^Q)/// = 137^^^(e2 — ei), and the QED first frequency of each 
spectroscopic series. 



k 


/. = e-' 


1373^2(62 - ei) 


wqed 


<1 


1 


185.99 


3.996x10-^ 


3.750x10-1 


7 


2 


307.63 


8.831x10^2 


6.944x10-2 


9 


3 


475.08 


2.398x10-2 


2.430x10-2 


11 


4 


577.99 


1.331x10-2 


1.125x10-2 


12 


5 


694.77 


7.667x10-=^ 


6.111x10-3 


13 


6 


826.22 


4.558x10-3 


3.685x10-3 


14 


7 


973.12 


2.790x10--'^ 


2.406x10-3 


15 


8 


1136.27 


1.752x10--'^ 


1.640x10-3 


16 


9 


1316.44 


1.127x10-3 


1.173x10-3 


17 


10 


1514.40 


7.403x10-^ 


8.678x10-^ 


18 


11 


1730.93 


4.958x10-^ 


6.600x10-^ 
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12 


1966.77 


3.379x10-^ 


5.136x10-^ 


20 


13 


2222.70 


2.341x10-^ 


4.076x10-^ 
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Caption to Table 1: Quantum number k of the QED transition k + 1 ^ k, numerically 
calculated angular momenta = in units of e'^/c, orbital frequencies in atomic units 

(137^)3 = 1373^2(^^2 — 61), the corresponding A;*^circular line of QED in atomic units, Wqed = 
Kp- — (^ipjp-) , and the integer q of Eqs. (45) and (47). 

Table 1 illustrates the fact that a resonance among PP modes predicts magnitudes pre- 
cisely in the atomic scale, as first discovered in Ref. [2]. In Ref. [2] we had to jump the 
integer q by twenty units for a complete quantitative and qualitative agreement with the 
Bohr atom. The qualitative agreement achieved by Table 1 is superior in this respect; but 
the correspondence is still not perfect for q < 7, after which q increases one by one, in qual- 
itative agreement with QED. The angular momenta in the second row of Table 1 should be 
compared with the QED values, i.e., = 137.036A;. In Table 2 we give the anomalous roots 
for q < 7 . Notice that the angular momenta are still of the order of Planck's constant, but 
the orbital frequencies do not correspond to any line of hydrogen. Since the approximation 
at Eq. (54) uses an expansion in 9, we should expect it to fail for the largest values of 9, as 
it did. 



20 





137^0^(69 - e^) 


a 


93.26 


3.16 


1 


51.06 


1.93 


2 


57.26 


1.36 


3 


76.14 


5.82 


4 


103.92 


2.29 


5 


140.37 


0.929 


6 



Caption to Table 2: Numerically calculated angular momenta h — ^ in units of e^/c, 
orbital frequencies in atomic units (137^)^ = 137^^^(e2 — ei), and the values of the integer q 
of Eqs. (45) and (47). 

In the following we give a second derivation of Poynting's resonance condition (52). The 
angular momentum of gyration along the orbital plane, calculated with Eq. (47) and aver- 
aged over the fast timescale is 

Ix + ily = /xr2|i?f ||i?f|^6iexp[i(e^-e? + ^)fit/^] (57) 



= ||i?f|exp[i(e^-e? + ^)Ot/^], 

9 

where we used Eqs. (1) and (2). The angular momentum vector of Eq. (57) rotates at the 
same frequency Awpp of Eq. (50), the ping-pong beat. Since angular momentum carries 
inertia, its rotation produces a gyroscopic torque on the orbit, so that the guiding-center 
motion should display an oscillation at the frequency Awpp of Eq. (50). The derivation of 
Eq. (57) assumed that the guiding-center orbit is a circular orbit of frequency Q, therefore 
we must have Q = Awpp. In Appendix C we give a third derivation of resonance (52), by 
averaging the electronic spin-radiation field over the fast timescale. The resonance is needed 
for the averaged spin-radiation field to rotate at the guiding-center frequency, so that it can 
participate in the guiding-center dynamics. 



VI. CONCLUSIONS AND DISCUSSION 

We stress that the QED condition, that 9 is of the order of the fine-structure constant, 
was never used in the calculations. The PP frequencies were calculated for arbitrary 9, and 
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it was the resonance condition of Poynting's theorem alone, a nonradiation condition, that 
calculated 9 in the atomic region. The quasi-degeneracy phenomenon was already found in 
Ref . [2] , and we expect it to be a universal property of electromagnetic dynamics with delay, 
as follows. The corrections to the limiting form (43) are controlled by powers of 1/A and 
9'^\, so that the resonance condition essentially poses a relation oc 1/A, which yields 
9 = c/|A| = c/(7rg), in agreement with Table 1, with Ref.[2] and with the Bohr atom of 
QED[3]. The finite- amplitude corrections discussed in Appendix C introduce again powers 
of p ~ l/l'^l; so that the qualitative behavior should persist. 

In Refs. [17] we argued that a stable orbit should emit at a neutrally stable frequency of 
its tangent dynamics. The physical process near the stable orbit is that the two terms on the 
right-hand side of Eq. (48) do not compensate exactly, so that there is a net radiation at the 
resonant frequency. Moreover, the guiding motion of the three-frequency orbit is not exactly 
a Coulombian circular orbit, because the perturbative scheme should correct the Coulombian 
approximation significantly, specially for the first states, where sup(p) ~ I/ttq' is larger. 
The orbital frequencies of Table 1 are only an approximation to the emitted frequencies, an 
approximation that should be worse for the low values of q, again in agreement with QED. 
More generally, according to Kurtzweil's small delays don't matter theorem [18], small 9 
orbits should be solutions of a limiting ordinary differential equation obtained by setting 
all delays to zero. This theorem yields that the limiting orbits should be the Coulombian 
orbits. In agreement with this, the gyration radius of Eq. (98) is a smaller fraction of the 
orbital radius for larger values of q , the sup(p) l/nq hmit of Appendix C. We call this 
limit the Kurtzweil correspondence limit. 

For 9 in the atomic magnitude, there is a nontrivial stabilization mechanism involving PP 
modes, a mechanism that selects discrete orbits by a resonance. The ping-pong modes form a 
channel that can interfere and absorb the energy radiated by the orbital motion, a dynamics 
that involves a resonance between the fast and the slow timescales. A multiscale solution, 
without assuming that the guiding-center motion is circular, is still to be worked out. The 
multiscale solution involves a general guiding-center slow motion plus a fast PP gyration of 
small amplitude. After the fast dynamics is balanced locally, the multiscale method should 
yield differential equations for the guiding-center orbit by a Fredholm alternativc[19-21]. 
Qualitatively, we expect the guiding-center equations of motion to accept richer orbits than 
illustrated in Figs. 1 and 3. The other possible orbital topologies are: (i) gyration plane 
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perpendicular to the normal connecting the particles, a dumbbell-like spinning, (ii) fast 
gyration restricted to the orbital plane, and (iii) non-planar guiding-center motion, like a p 
orbital. The intuition about the dynamics should follow the general guidelines learnt from 
analyzing the circular-guiding-center-three-frequency-orbits of Figs. 1 and 3, as summarized 
in the following. For a resonant orbit, the electronic spin-radiation force, averaged over the 
fast timescale, yields an 0(1) slowly rotating force at the proton that perturbs the protonic 
equilibrium position along the guiding-center orbit. The protonic attraction cancels the 
electronic self-interaction, so that the guiding-center motion does not dissipate (any guiding- 
center motion). The averaging of the electronic spin-radiation force involves a resonance, 
and this is the reason for discrete orbits from the perspective of the equations of motion for 
the guiding-center. We derived the resonance condition in three different ways here; (i) In 
Section V we derived the resonance condition using Poynting's theorem, (ii) Still in Section 
V, we derived the same resonance condition using general arguments about the gyroscopic 
torque, and last (iii) In Appendix C we derived the resonance condition from the detailed 
balancing of the guiding-center dynamics. 

In the process of emission, the sharp line is emitted when the dynamics is approaching 
the stable orbit, as discussed in Rcf. [17]. This long timescale is to be compared with 
the time of spontaneous decay of QED for hydrogen; about 10^ orbital periods or 10~^° 
seconds. On the other hand, the dynamics has a fast timescale, the frequency of the PP 
modes, of about nq/O ~ 1000 times the orbital frequency. This fast frequency is about 
10^° Hertz and resonates with the X-ray frequencies used in the Compton effect [22], an 
interesting coincidence. The ping-pong frequencies of hydrogen have the magnitude of the 
zitterbewegung frequency of Dirac's relativistic version of Schroedinger's equation[24]. In 
modern Aharonov-Bohm experiments [23], the ballistic electron passes at a manometric 
distance from the electrons inside the solenoid, such that the PP frequencies, that fall as 
1/r, can be even at X-rays. The magnetic field interacts with the PP oscillations of the 
solenoid electrons, which in turn influences the baUistic electron. This fact that electrons 
can play ping-pong-at-a-distance with frequencies up to X-ray suggests the need of an X-ray 
shielding for the solenoid of the Aharonov-Bohm experiments, e.g., a layer of lead. 

The three-frequency orbit solves several conundrums of the hydrogen atom, paradoxes 
that were created by guessing that the equations of classical electrodynamics would accept 
non-stiff planetary-like orbits at every scale. The ping-pong phenomenon is a non-trivial 
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feature that is not present in the ordinary differential equations of planetary dynamics. The 
qualitative agreements with QED are listed in the following; (i) the resonant orbits are 
quantized by integers that appear naturally because of the delay, and the orbital frequencies 
agree reasonably with the Bohr circular lines, (ii) the angular momenta of the resonant 
orbits turn out to be approximate integer multiples of Planck's constant, (iii) The emitted 
frequencies are given by a difference of two eigenvalues of a linear operator, the frequencies 
of the PP modes, analogously to the Rydberg-Ritz combinatorial principle of QED. (iv) 
the averaged angular momentum of gyration is of the order of the electronic spin- angular- 
momentum of QED [26, 27]. The approximation for the gyration angular momentum is not 
satisfactory yet though. 

Wc exhibited a new orbit of the hydrogen atom of Dirac's electrodynamics oi point charges 
with retarded-only fields. The Lyapunov stability of this orbit should be further investigated. 
Because we are dealing with a physical theory, the equations should be sufficiently robust 
to allow some kind of numerical integration [28] . The stability of the three- frequency orbit 
poses a linear set of delay equations, a dynamical system that needs an initial function as 
the initial condition, just like Schroedinger's equation. It would be interesting to learn if this 
linear operator has a self-adjoint Fredholm alternative [19, 20], like Schroedinger's equation 
does. The frequency of PP modes is proportional to the inverse of the interparticle distance, 
which turns out to be proportional to the electrostatic potential energy (an accidental anal- 
ogy, because these are different phenomena). The potential energy is thought to be the 
reason why the 1/r appears in Schroedinger's equation; and it is interesting to notice that a 
linearized equation for PP modes has the same generic form. The ping-pong modes appear 
universally in electromagnetic many-body dynamics because of the delay, a phenomenon 
that has been so far overlooked. The interference mechanism of Poynting's theorem can 
produce orbits that do not radiate, if a resonance is satisfied. This resonance turns out to 
be satisfied precisely in the atomic magnitude, a surprise that makes this dynamics inter- 
esting for theoretical physics. We hope that our preliminary findings serve to motivate and 
guide further studies of this dynamics, and there is much to be settled yet. 
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VIII. APPENDIX A: ELECTRODYNAMICS OF POINT CHARGES 

In 1947 Eliezer generalized Dirac's covariant subtraction of infinities [10]. In this general- 
ized electrodynamics of point charges [10], the field produced by the point charge is supposed 
to be the retarded field plus an intrinsic free field G 

f; = i^;,.et + Gi- (58) 

The free-field G used by Ehezer in Ref. [10] is finite along the particle's trajectory and 
vanishes when the particle is at rest, 

G; = r(F;,.et-^;,aJ, (59) 

where F is a constant [10]. This generalized electromagnetic setting is henceforth called the 
Eliezer setting (ES). In Eliezer's theory, the electron and the proton of a hydrogen atom, of 
charges —1 and -|-1 respectively, have the following equations of motion [10] 

m^v,, = ^(1 + 2F)[i>^i - \\vA\^v,,] - [F;i,,„ + (1 + F)F;2,,et - ^F^.^adv]^,,, (60) 
m2V,2 = + 2T)[v,, - \\v2\\'v,,] + [F;^,, + (1 + F)^;,^,,, - FF;,^„,Jz;,2, 

where the v^i stand for the four components (// = 1,2,3,4) of the Minkowski velocity of 
particle i, double bars stand for the Minkowski norm and the dot indicates derivative respect 
to the proper time of each particle. Equation (60) includes an external field -F^j,„ produced by 
the other charges of the universe at each particle's trajectory, which vanishes for the isolated 
electromagnetic two-body problem. The first term on the right-hand-side of Eq. (60) is the 
sourceless combination of half of the retarded Lienard-Wiechert self-potential minus half 
of the advanced Lienard-Wiechert self-potential, all multiplied by (1 -|- 2F) and evaluated 
at the particle itself [10]. This term generahzes Dirac's self-interaction force[l]. The last 
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two terms on the right-hand side of the first fine of Eq. (60) define the interaction with 
the retarded potentials of the proton at the electronic position, F^2,refi interaction 
with the advanced potentials of the proton at the electronic position, Fj^^adv 
four interesting limits; (i) F = is Dirac's electrodynamics with retarded-only fields [1], 
(ii) r = — 1 is a non-physical theory with advanced-only interactions, (iii) F = —1/2 is the 
action-at-a-distance electrodynamics [29] ( notice that the self- interaction vanishes), and 
(iv) The limit when F tends to —1/2 from above is the dissipative Fokker setting of Ref. [2], 
with a charge renormahzation controlled by (1 -|- 2F). 

For the non-specialist reader, in the following we write the equations of motion of Dirac's 
electrodynamics, F = 0, in the intuitive form of physics textbooks. For the isolated hydrogen 
atom, the spatial component of the electronic equation of motion, the first line of Eq. (60), 
multiplied by ^/l — |viP, yields 

m,^{-=Z^) = Fi - E2 - X B2, (61) 
ah y/l - |vi|2 

where Vi is the Cartesian electronic velocity[12, 25]. The equation for the proton is obtained 
exchanging the indices in Eq. (61) and multiplying the last two terms on the right-hand 
side of Eq. (61) by minus one, to account for the positive protonic charge. The first term on 
the right-hand side of Eq. (61) is the force Fi of the electronic fields on the point electron 
itself, and is called the Lorentz-Dirac self-force, 

Fi = '^llisii + 7?(vi • ai)vi + 3-fl{Yi ■ ai)[ai + ^l{vi ■ ai)vi]}, (62) 

where 7^ = 1/ -^1 — |vip and ai and ai stand for the Cartesian electronic acceleration and 
time-derivative of the Cartesian electronic acceleration, respectively (cf. page 116 of Ref 
[25]). In Eq. (62) a dot between two Cartesian vectors indicates scalar product and a dot 
over a vector indicates derivative respect to laboratory time ti. The second term on the 
right-hand side of Eq. (61) is the electric force, i.e., the electronic charge, e = —1, times 
the electric field of the proton acting on the electron, E2. In our unit system with c — 1, 
the retarded electric field of the proton, of charge +1, is given by the Lienard-Wiechert 
formula[15] 

J, ^ ni26 - V26 _^ ni26 X [(ni26 - V26) X 3i2b] /ggv 

^ llbi^ - ^12b-y2bTrl2b (1 - ni26-V26)3ri26 

where V26 and a26 are respectively the Cartesian velocity and Cartesian acceleration of the 
proton at the retarded time t2b and ^2b = I/a/I — |v26P. In Eq. (63), unit vector n^b 
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connects the retarded protonic position to the present electronic position, and r^b is the 
interparticle distance along the retarded light-cone. The advanced fields are obtained by 
replacing by — V26 in Eq. (63), and are not used here, since we deal we Dirac's retarded- 
only theory. The third term on the right-hand side of Eq. (61) is the magnetic force, i.e., the 
electronic charge times the vector product of the electronic velocity Vi by the magnetic field 
B2 of the proton. The retarded magnetic field of the proton is given by the Lienard-Wiechert 
formula [15] 

B2 = ni26 X E2. (64) 

For the advanced magnetic field, formula (64) includes a minus sign, so that the Poynting 
fiux E2 X B2 calculated with the purely advanced fields of an unperturbed circular orbit is 
an incoming energy flux. The combination — (E2 -|- Vi x B2) is usually called the Lorentz 
force, not to be confused with the Lorentz-Dirac self-force (62). The electric and magnetic 
fields (63) and (64) depend on the retarded position, velocity and acceleration of the proton, 
and this is where delay enters in the dynamics. Equation (62) can be rearranged as 

Fi = l^^U^{k^ + 37?(vi • ai)ai}, (65) 

where we introduced the dyadic matrix Ui = /-|-7^ViV^, a non-singular matrix whose inverse 
is U^^ — I — ViV*. The left-hand side of Eq. (61), i.e., the variation of the momentum, can 
be expressed using Ui as 

^^i^Vi-ViP^^'"^'^^^^^" ^^^^ 

Using Eqs. (65) and (66), equation of motion (61) can be expressed as 

miai = ^7i{ai + 37^(vi ■ ai)ai} - ^/l^^J^U^^E^ + vi x B2), (67) 

where E2 and B2 are the electric and magnetic fields of the proton, respectively, and we 
have cancelled the invertible matrix Ui and a power of 7^. Equation (67) has a familiar 
Newtonian-likc form, but it is still the full rclativistic equation. In the low-velocity limit, 
the first term of the right-hand side of Eq. (67) reduces to the third derivative of the position 
multiplied by 2/3, which is called the Abraham-Lorentz-Dirac self-force. 

There exists an action formahsm for the Lorentz- force sector of Eq. (61), i.e., Eq. (61) 
without the first term on the right-hand side. We give a general action including the pa- 
rameter r solely for the calculations and crosschecking of Sections III and IV. A reader 
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interested only in Dirac's theory can set F = 0. This action involves only the instantaneous 
position and velocity of particle 1, and is composed of a kinetic term plus the following 
electromagnetic action, 

e = -r / - • dt, + (1 + r) / " ' "'^^ . dh. (68) 

J ri2a{l + ni2a " V2a) J ri2b{l - U^b ■ ^2h) 

In Eq. (68), Vi stands for the Cartesian velocity of particle 1 at time ti, while and V2b 
stand for the Cartesian velocities of particle 2 respectively at the advanced time t2a and at 
the retarded time ^26- The vector ni2a in Eq. (68) is a unit vector connecting the advanced 
position of particle 2 to the position of particle 1 at its present time while unit vector 
ni26 connects the retarded position of particle 2 to the present position of particle 1 at time 
ti. Still in Eq. (68), ri2a and ri26 indicate the interparticle distance along the advanced 
and retarded light-cones, respectively. To derive the equations of motion of particle 1, one 
needs to add the usual kinetic action to Eq. (68), Ki = — J m\^Jl — \'Vi\'^dti ( the integral 
of the kinetic Lagrangian). In Ref.[12] it is shown that formal minimization of the sum of 
action (68) with the kinetic action yields the equations of motion of particle 1 suffering the 
electromagnetic fields of particle 2, i.e., Eq. (61) without the self-interaction term. Each 
integrand of the right-hand side of Eq. (68) is a familiar electromagnetic Lagrangian 

Lc ^ ^^7^'7":\ ^ -{Ve - vi • A,), (69) 

where Vc and Ac are the Lienard-Wiechert scalar potential and the Lienard-Wiechert vector 
potential, respectively. We introduced the quantity c = ±1 in Eq. (69) to indicate if the 
interaction is along the advanced or the retarded light-cone. Equation (69) with c — 1 
indicates that particle 2 is in the future, while c = — 1 indicates that particle 2 is in the 
past. The quantities of particle 2 in Eq. (69) are to be evaluated at the time t2c defined 
implicitly by 

t2c = + — , (70) 

c 

where ri2+ is the distance along the advanced light-cone and ri2- is the distance along the 
retarded light-cone. According to Eq. (70), the time lag along the advanced light-cone is 
ri2+ and the time lag along the retarded light-cone is — ri2_. 

The shortest way to obtain the equations of motion of Dirac's electrodynamics in any 
given coordinate system is to start from the Euler-Lagrange equations of action (68) plus the 
relativistic kinetic action. This yields the dynamics without self-interaction of each particle 
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suffering the electromagnetic fields of the other particle. The self-force can be added to the 
Euler- Lagrange equations, watching carefully for the correct multiphcative factor. The stiff 
limit of PP modes is determined by the largest-order derivative appearing in the linearized 
equations of motion. In this limit, the contribution of the self-interaction to the linearized 
dynamics about a circular orbit is given by the Abraham-Lorentz-Dirac self-force 

Frad = + 2T)k. (71) 

The contribution of the other nonlinear terms, that become important at a finite distance 
from the circular orbit, is discussed in Appendix C. The electrodynamics of point charges 
is discussed at length in Refs. [10, 25], while the Lienard-Wicchert potentials, Lienard- 
Wiechcrt fields, Poynting 's theorem and the physics of electrodynamics is found in innu- 
merous textbooks, e.g. Refs. [12, 15]. 



IX. APPENDIX B: VARIATIONAL EQUATIONS ALONG THE z DIRECTION 

In this appendix we derive the linearized variational equations for displacements perpen- 
dicular to the orbital plane, henceforth called the ^-direction. Since Zk = is an exact 
solution of the equations, in this Section we are doing Lyapunov stability. The variational 
dynamics along the z-direction is uncoupled from the planar variational equations up to 
the linear order. In the same way of Section IV, we expand to second order the implicit 
hght-cone condition and action (68). The Cartesian coordinates of a transversely perturbed 
circular orbit are defined as 

Xk + iyk = rbdkexp{iQt), (72) 
Xk-iyk = rbdlexp{-iQt), 
Zk = nZk, 

where k = 1 for the electron and k = 2 for the proton, is the transverse perturbation, 
di = bi and d2 = — &2 are defined in Eq. (4) and fl is the orbital frequency (1). In the 
following we calculate the delay function (f^ of (13) by expanding the light-cone time t2c 
about the constant lag r?, up to the second order in Zi and Z2. The distance r^c in Eq. 
(70) is evaluated from the position of particle 1 at time ti to the position (72) of particle 2 at 
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time t2c- Using t2c defined by Eq. (70) and orbit (72), this implicit distance r^c — |^2c — ti\ 
evaluates to 

rL = rlil + (t>f = rl[hl + 6^ + 26162 cos((^ + Oc)] + rl{Z^ - Z^^f, (73) 

where we expressed </7p in terms of the scaled function of (24). The Z variations decuple 
from the planar variations because powers of Z always appear squared, so that there is no 
mixed linear term of Z times a linear perturbation of the planar coordinate in Eq. (73). 
Expanding Eq. (73) up to the second order and rearranging yields 

02 + 250 = (Zi - Z^cf. (74) 

Equation (74) is a quadratic equation for 0, and the regular solution is given, up to second 
order, by 

0=^(Zi-Z2e)2. (75) 

The coordinate Z2 appears evaluated at the advanced/retarded time in Eq. (75), and to 
obtain the action up to quadratic terms it is sufficient to keep the first term Z^ijx + + </?) 
~ Z^ijx + cB) = Z^c- Using orbit (72) to calculate the numerator of Lagrangian (69) yields 

(1 - Vi • V2c) = 1 + C0S(^ + 0(^)6162 - ^'^1^2c ~ (76) 

c - e\s - 1)0 - e^ZiZ2c, 

while the denominator of Lagrangian (69) evaluates to 

ri2c(l + ^^1?^^) = n[l + + echb2 Sm(ec + (^J + ec(Z^ - Z2e)^2c]- (77) 

c 

Notice that the quadratic term Z2cZ2c on the right-hand side of Eq. (77) can be dropped 
because it represents an exact Gauge that does not affect the Euler-Lagrange equations. We 
also expand the argument of the sine function of the right-hand side of Eq. (77) until the 
linear term in 93^, so that the quadratic approximation to Eq. (77) is 

ri2c(l + '"'''""''' ) ~ n\S + C0 + ^cZiZ2e], (78) 
c 

where the equivalence sign henceforth means equivalent up to a Gauge term of second 
order. Even if a quadratic Gauge term appears in the denominator, in an expansion up to 
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quadratic order it would still produce a Gauge and therefore it can be dropped directly from 
the denominator. The Lagrangian of action (68) expanded up to second order is 

£t ~ (-^){(i + r)[i - e^cs^z.z^- - ^{z, - z^-f + ec^sz.z^^] (79) 
-r[i - d^cs^z^Z2+ - ^(Zi - Z2+f - ec^sziZ2+]}. 

We henceforth disregard 0(9'^) corrections, and substitute C = 5* = 1 in the coefficients of 
Lagrangian (79). Last, we need the kinetic Lagrangian along orbit (72), i.e.. 



Ti = -mi - v! = -^Vl - ll&'Zl (80) 



where the dot means derivative with respect to the scaled time r, 7^ ^ = \/l — vf , and we 
have used that flrb = 9. Equation (80) expanded up to the second order is 

T. = (i){==p + |z? + ...}, (81) 

where ei = r^mi7iQ^. Using Eqs. (1) and (2), we obtain ei = mi/ii — M/m2. The equation 
of motion for particle 1, without the self- interaction term, is the Euler-Lagrange equation 
of the quadratic Lagrangian 

L^y^ = Ti + £r. (82) 
The Lagrangian sector of the equation for Zi is 

eiZi = -[Zi - (1 + r)Z2- + rz2+] + 9[rz2+ + (i + r)Z2_) + 9^[rz2+ - (i + r)Z2_]. (83) 

Notice that the left-hand side of Eq. (83) can be written as 

e,Z, = rlm,^,n'Z, = rf-^, (84) 

which is the force along the z-direction multiphed by r^. As explained above Eq. (40), to 
account for self-interaction we must add to Eq. (83) the Abraham- Lorentz-Dirac self-force 
(71) multiplied by , 

rl¥rad = \{l + 2T)9^Zi. (85) 
The full linearized variational equation for Z^ is 

eiZi = ^(1 + 2r)^3^i - [Zi + rz2+ - (1 + r)Z2_] (86) 
+9[TZ2+ + (1 + r)Z2_] + ^'[rz2+ - (1 + r)Z2_]. 
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The linearized equation for Zi is obtained by interchanging Z\ by Z2 and e\ by 62 in Eq. 
(86). Comparing Eq. (86) to Eq. (30) of Ref. [2] we find that Eqs. (29) and (30) of Ref.[2] 
are both missing a factor in front of the Z\ term, which is a typo. After Eq. (30), the 
other equations of Ref. [2] have the self-interaction included correctly. 

A ping-pong normal mode is obtained by substituting Zi = Aexp(A2f2t/^^) and Z^ = 
B exp{Xzflt / 9) into Eq. (86) and the corresponding linearized equation for the proton. 
Again we use a general A^, but a harmonic solution needs a purely imaginary A^. Setting 
the determinant to zero yields 



i + ^-i(i + 2r)A^ G{e,K 



mie^ 3 

X ^ 1 -U 



G{B.K) l + S^-|(l + 2r)A^ 



0, (87) 



where G{Q, A^) = [1 - (1 + 2r)A^ - A^] cosh(A^) - [(1 + 21) - A^ - (1 + 2r)A3 sinh(A^). The 
stiff-limit is when |A^| is large, and we should keep in mind that the hyperbolic functions in 
G{6,Xz) can acquire a large magnitude [11]. Multiplying determinant (87) by iJ,6'^ / {M X'^) 
we obtain 

1-1(1 + 2r)0'X, + i A9*A^ (88) 



z 

i2n 



up to small 0{d ) terms. The stiff-mode condition defined by Eq. (88) with V — —1/2 is 
equation (33) of Ref. [2], i.e.. 



1 - + ^Q'^l - ^[(1 - ^)[cosh^(A.) - ^sinh(2A.)]^ = 0. (89) 

In Ref [2] there is another typo in passing from Eq. (33) to Eq. (34); Eq. (34) is missing 
a bracket that should start after the {pB"^ jM) factor and close at the end of Eq. (34). The 
stiff-hmit in Dirac's theory with retarded-only fields, (F = 0), is 

and the appearance of the negative exponential is related to the retardation-only. Comparing 
Eq. (88) to Eq. (41) we find that the quasi-degeneracy phenomenon exists only for (i) F = 
, i.e., Dirac's theory, (ii) F — —1, a non-physical advance-only case, and (iii) F = — 1/2 , the 
action-at-a-distance electrodynamics[29] and the dissipative Fokker dynamics of Ref. [2]. 
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X. APPENDIX C: THE PING-PONG SOLUTIONS 



In this Appendix we discuss the existence of a harmonic orbit such as given by Eq. (47) 
and illustrated in Figs. 1 and 3. In the following we show that an orbit such as (47) exists 
in the limit where p^. = sup(|i?^|, l-R^^I) ^ j^, a limit where the stiff gyration approaches 
the speed of light. In this hmit, the dominant field of the particle is the far-electric field, 
Eq. (48). The near-electric field is the first term on the right-hand side of Eq. (63), which 
vanishes because oo. The far-magnetic interaction is the next-to-leading term in size, 

and is disregarded in the following. We also disregard the self-interaction, Eq. (62), because 
its contribution to the force normal to the velocity is smaller. The equation of motion for 
the proton is obtained by exchanging indices in Eq. (67), and the leading fast dynamics 
normal to the fast velocity is 

(1 - n • Vi_)3r2i6 

where underscore minus indicates that particle 1 is in the past light-cone. The unit vector 
n points from the retarded position of particle 1, which is not indicated with underscore to 
avoid an overloaded notation. For the PP modes, the fact that the electronic coordinates 
on the right-hand side of Eq. (91) are evaluated at ti- ~ t2 — makes a lot of difference, 
because the PP modes execute complete periods during this time-lag. Since the gyration 
amplitude is small compared with r;,, we henceforth assume that the distance r^ib along the 
retarded light-cone on the right-hand side of Eq. (91) is constant and given by r^. Along 
the orbit of Fig. 1, the fast gyration can be parallel to the normal, so that (1 — n ■ vi_)^ 
becomes arbitrarily small, and we henceforth approximate the denominator of the right- 
hand side of Eq. (91) by (1 — |A|pi)^r6. We also approximate the square-root on the 
left-hand side of Eq. (91) by (1 — \\\^p'iY^'^- The stiff-limit is obtained by substituting orbit 
(47) with py. — sup(|ii!^|, l-R^^I) into Eq. (91), and taking the Fourier component of Eq. 
(91) along the PP frequency wpp ~ 'KqVt/d. We also multiply Eq. (91) by r^, so that the 
order-of-magnitude balancing posed by the Fourier-transformed version of Eq. (91) is 

(l-|A|Vi)^/ ^ 

^,e'P'"" {i-\\\p,Y 

The equation of motion for the electron is obtained by interchanging indices in Eq. (92). A 
salient feature of Eq. (92) is that A must be purely imaginary and multiple of tt, A = T:qi, 
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2^2 - — VTTT^Pi exp(-A). (92) 



so that the p's can be real in Eq. (92), i.e., the phase-shift 2Im(A) must be a multiple of 
27r. Considering the fraction on the right-hand-side of Eq. (92) as a given number, Eq. (92) 
and the corresponding equation for the electron are two linear homogeneous equations for 
and Equating the determinant of this homogeneous system to zero and disregarding 
smaller terms yields 

/.gSxp(-2A) ^ (1 - \\\p,f'\l - \\\p,fl' 3 
M V(l + |A|Pi)(l + |Ah) ■ 

Equation (93) gives the main finite-amplitude correction to Eq. (43), to which it reduces 
for Pi = P2 = 0. The finite-p corrections cancel the a's, which was the justification for the 
approximation explained above Eq.(54). When sup(p]^,p2) approaches from below, the 
right-hand side of Eq. (93) matches the O{0'^) left-hand side, so that Eq. (93) accepts a 
purely harmonic solution with A = Trgi, as we wanted to demonstrate here. 

After solving for the ping-pong oscillation, the dynamics at the slow frequency is the 
next-to-leading of the multiscale solution. For the electronic motion, Eq. (67) has a zero- 
order dissipativc force along the unperturbed orbit, as illustrated in Fig. 3. This offending 
force has a component along the electronic velocity that is given by 

rlFr = rf-x, = -%^'rlx, = (94) 

where we multiplied the force by r^, as in Eqs. (40) and (85). The Lorentz- force of the 
proton along the unperturbed circular orbit is almost normal to the electronic velocity and 
does not contribute much for the dissipation. On the equation of the proton, the main 
offending force at zero-order is the delayed interaction with the electron, instead of the 
much smaller protonic self-interaction. Using Page's series in the same way of Ref. [17], we 
find an offensive zero-order force against the protonic velocity of the same magnitude of force 
(94). It turns out that it is impossible to cancel the zero-order force (94) with the linear 
terms of the variational equations, which is shown by averaging the linearized equations over 
the orbital period, an averaging that yields a 4 x 4 linear system with no solutions. The 
balancing of the offensive zero-order force is nevertheless possible by using the resonance 
between ping-pong modes that is introduced at quadratic order by the spin-radiation term, 
as discussed below. 

Multiplying the offensive force (94) by the velocity along the circular orbit yields an 0{9^) 
dissipated power. Again, the discussion of the dissipated power must start from the leading 
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fast dynamics, as follows. Taking the scalar product of Eq. (67) with Viand disregarding 
the contribution of the protonic fields for the dissipated power yields 

2 

mi(vi • ai) = -7i{(vi • ai) + 37i(vi • ai)^}. (95) 

The bracket on the right-hand side of Eq. (95) is multiplied by the possibly large factor 7^^, 
so that it must vanish in the large-|vi| limit, i.e., 

(vi-ai)2^--^(vi-a0. (96) 

The offending force (94), multiphed by the slow guiding-center velocity 9 gives only an 0{9'^) 
contribution to the right-hand side of Eq. (95). The dominant contribution is given by the 
PP oscillation, of frequency fi|A|/^, and defined by Eqs. (46) and (47). We henceforth 
replace the ai on the right-hand-side of Eq. (96) by — (r2^|Ap/^^^)vi. Multiplying Eq. (96) 
by rl and using the fast component of orbit (47) with = sup(|i?f |, l-Rf^l) yields 

|A|V?cos2(x)c.i(l-p?|A|2), (97) 

where x is the angle between the velocity and the acceleration. In the limit when |vi| 1, 
this angle x must be close to ninety degrees, since |vi| must be less than one and (vi • ai) = 
^(■'^y-). Using Eq. (97) together with | cos(x)| < 1 yields a lower-bound for 

Next-to- leading in the multiscale solution is the guiding-center dynamics, i.e., the equilibra- 
tion of ofTending-force (94), a balance that takes place after the fast dynamics is established 
as a harmonic oscillation with a radius given by Eq. (98). In the following wc give a third 
derivation of resonance (52) by averaging the equations of motion over the fast timescale. 
Due to the larger protonic mass, we henceforth assume that the protonic field at the elec- 
tron averages essentially to the Coulombian field, a force that rotates at the guiding-center 
frequency. On the other hand, in the limit Pi\X\ 1, the electronic field on the proton, 
averaged over the fast timescale, is significantly changed. Analogously to the derivation of 
Eq. (57), the electronic spin-radiation field, second term on the right-hand side of Eq. (91), 
averaged over the fast timescale, rotates at the guiding-center frequency fl if resonance (52) 
holds. This so averaged field adds to the rotating electronic Coulombian field (not included 
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in Eq. (91)), so that the averaged attractive force on the proton deviates from the diametral 
direction. The proton repositions along the circular orbit, until its velocity becomes per- 
pendicular to this perturbed centripetal force, as illustrated in Fig. 3. After repositioning, 
the angular distance between particles at the same time is less than 180 degrees. The self- 
consistent repositioning stops when the Coulombian field of the proton along the electronic 
velocity acquires a component to balance the electronic self-interaction (94). The perturbed 
equilibrium position along the guiding-center orbit is illustrated in Fig. 3. The situation 
of Fig. 3 is possible only at resonance, when the averaged spin-radiation force acting on 
the proton rotates at the guiding-center frequency. In Ref. [17] we also used a resonance 
between mutually orthogonal vibration modes. The present work goes beyond our simple 
estimates of the helium dynamics of Ref. [17]. Unlike the Coulombian orbits of helium[17], 
Coulombian circular orbits of hydrogen radiate in dipole, and we have seen here that the 
fast PP oscillations are essential to cancel this dipolar radiation, a complex dynamics that 
demands a multiscale solution. The PP modes are lost when the delay is expanded, as in 
Page's series of Ref. [17], a nontrivial manifestation of stiffness. 

For 9 in the atomic magnitude, lower-bound (98) is a few percent of the interparticle 
distance, nevertheless the distance p^r^ is already some 1000 classical electronic radii, much 
larger than the radius of the fat Lorentz electron[27]. To improve on estimate (98) one needs 
the relation between cos(x) and (1 — p^|A|^), which is beyond the present work. To solve 
Eq. (93) it suffices that either ]i?f] or \R^\ approaches the upper limiting value the 
other Rk can take a much lower value, i.e., a solution exists already in the limit when just 
the electron rotates near the speed of light. As an estimate for the angular momentum of 
gyration, Eq. (57), we use \R\\ and \R\^\ given by Eq. (98), yielding 

47r& q 

Calculating angular momentum (99) with the first hue of Table 1 gives \lxy\ — 213.5. The 
angular momenta estimated by Eq. (99) depend on g, but are in the correct order of 
magnitude. The electronic spin-angular-momentum of QED is independent of the orbital 
quantum number and given by \s\ = \/3^/2 = 118. 
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XI. CAPTIONS 



Fig. 1: Guiding-center circular orbit (dark lines) and particle trajectories gyrating about 
the guiding-center (grey lines) for a three-frequency orbit near resonance. Illustrated is also 
the beat of the electron's gyration radius at about the orbital frequency. Drawing not on 
scale; the beat of the protonic gyration radius is not illustrated. Arbitrary units. 

Fig. 2: Coulombian circular orbit with particles in diametral opposition at the same time 
of the inertial frame (empty circles). Also indicated is the retarded position of both particles 
(solid circles), and the angle travelled during the light-cone time. Drawing not on scale; The 
circular orbit of the proton and the retardation angle have an exaggerated magnitude for 
illustrative purposes. Arbitrary units. 

Fig. 3: Perturbed guiding-center orbit, with particles no longer in diametral opposition 
at the same time. Positions at the same time are the solid circles, the angular distance at the 
same time is (tt — 5). At the retarded position of the proton (empty circle), averaged force 
Fi2 is perpendicular to the guiding-center velocity of the proton. At resonance, averaged 
force Fi2 rotates at the guiding-center frequency. Protonic orbit and self-force have an 
exaggerated magnitude for illustrative purposes. Arbitrary units. 
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